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Abstract 

■ Recently Murayama and Nomura proposed a simple scheme to construct the 
gauge mediation models, using metastable supersymmetry breaking vacua. It has 

C^l a possibility to predict the ultra light gravitino mass < 16 eV, while such a 

(— I ' light gravitino may destabilize the metastable vacua. We investigate stability of the 

metastable vacuum of their model. The transition rate from the false vacuum to 
true ones is evaluated by numerical calculation, including the Coleman- Weinberg 
J> ' potential destabilizing the metastable vacuum. It is found that when the messenger 

■ sector is minimal, stability of the metastable vacuum imposes an upperbound on 
^f) '. squark mass Mg for the ultra light gravitino as Mg < 1800 GeV at most. Squarks 

I with this mass may be found in the LHC experiments. 

od 

o 

o 

> 

X 



1 Introduction 



Introduction of Supersymmetry (SUSY) is one of the most attractive directions to con- 
struct models beyond the Standard Model (SM), since SUSY cures the hierarchy problem. 
However, when SUSY breaking mass terms of squarks and sleptons are generic in the 
SUSY SM, too large flavor changing neutral current (FCNC) processes are predicted. 
This problem is called as the SUSY flavor problem. In addition, the SUSY SM also has 
the SUSY CP problem. The SUSY breaking trilinear and bilinear scalar coupling terms 
{A and B terms) and the gaugino mass terms have CP phases generically. The predictions 
for electric dipole moments are rather large, though the constraints by experiments are 
severe. 

The Gauge Mediated SUSY Breaking (GMSB) models [H |2] are an excellent solution 
to these problems. SUSY breaking in a hidden sector is transmitted to the visible sector 
through the SM gauge interactions, which are flavor-blind. The SUSY breaking mass 
terms of squarks and sleptons are then flavor-blind, so that dangerous FCNC processes 
are naturally suppressed. In addition, the phases of the A and B terms are automatically 
aligned to those of the gaugino mass terms in the Minimal Gauge Mediation (MGM) model 
[3j, since the A and B terms are vanishing at tree level in the model. In the case, the 
SUSY CP problem is avoided^^ . 

In the GMSB models, the lightest SUSY particle is the gravitino, with mass —(1 
- 10 GeV). An ultra light gravitino 777.3/2 ^ 16 eV is known to have no cosmological grav- 
itino problem, such as overclosure [5] and too large free-streaming lengths [6]. However, 
only a few models were previously proposed in order to realize the ultra light gravitino [7]. 

Recently Murayama and Nomura proposed in Ref. [8] a simple scheme to construct the 
GMSB models along a metastable SUSY breaking scenario. Metastable SUSY breaking is 
recognized as an interesting alternative SUSY breaking scenario, after Intriligator, Seiberg 
and Shih found that even simple vector- like theories, such as SUSY SU(N) QCD, can break 
SUSY on metastable local minima [9]. Many phenomenological models with metastable 
SUSY breaking vacua are nowadays discussed. 

The scheme proposed by Murayama and Nomura has a quite simple structure and a 
possibility of the ultra light gravitino mass 7773/2 ^ 16 eV. However, the stability of the 
metastable vacua becomes weaker for smaller gravitino mass. It is not clear whether the 
lifetimes of the metastable vacua are longer than the age of the universe. 

In this paper, we study stability of the SUSY breaking metastable vacuum in the 
GMSB model which is realized in their scheme as an effective theory jSj , in the case of the 
ultra light gravitino mass. The lifetime of the metastable vacuum is numerically evaluated 
in the semiclassical techniques [10]. In the evaluation we include the Coleman- Weinberg 
(CW) potential in the scalar potential, which destabilizes the metastable vacuum. In 
Ref. [8], the authors estimated the lifetime by a triangular approximation [11]. However, 

It is shown in Ref. [4J that the latest experimental values of the anomalous magnetic moment of the 
muon and the branching ratio oi B Xs'y favor a low messenger scale 100 TeV) in the MGM model. 
This is quite consistent with the ultra light gravitino discussed in this paper. 
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this approximation is not accurate especially in this model, since several scalar fields are 
relevant to the phase transition, as the authors mentioned in their text. In the case, 
numerical calculation is required to estimate a reliable value of the transition rate. 

We also derive upperbounds on the squark mass and the cutoff scale for ultraviolet 
(UV) models for ^ 16 eV. The bound on the gravitino mass can be translated 

to those on the squark and gluino masses in this model, since the visible sector is more 
directly coupled with the hidden sector in this model. 

This paper is organized as follows. In the next section, we review the model discussed 
in Ref. [8]. In Section 3, we discuss the stability of the metastable vacuum and constraints 
on parameters of this model. Section 4 is devoted to summary. 

2 Model 

In this section we review the GMSB model with the metastable vacuum discussed in 
Ref. [8]. This model is realized in their scheme to construct the GMSB model, and has a 
quite simple structure. The relation between the gravitino mass and the stability of the 
metastable SUSY breaking vacuum is also clarified here. 

The superpotential and Kahler potential in this model are given by 

H/ = _^2^ + ^5//_M//, (1) 

^=i/r+i/r+i^r-^+c?(^), (2) 

where 5* is a gauge singlet chiral field and / and / are messenger fields, which transform 
as 5 and 5* under the SU(5)gut gauge symmetry, respectively. Here, /i, k and M can be 
taken to be real and positive without loss of generality, which are taken as free parameters 
in this model. 

This model is an effective theory below the UV cutoff scale A. It is assumed that higher 
terms of 15*1 appear after some particles with cutoff-scale mass are integrated out. The 
negative sign of the \S\'^ term is required for deriving a local minimum near the origin. The 
expansion in the above Kahler potential implies that 5* has some approximate conserved 
charge. For example, we can assign f/(l)_R charges of 2, and to S*, / and /, respectively, 
although the last term of Eq. ([1]) violates this symmetry explicitly. 

The tree-level scalar potential is given from Eqs. ([1]) and ([2]) as 

+ \KS-M\\\f\' + m + Y.\9^Dj. (3) 

i 

Here, we use the same notation for scalar fields as superfields, and Di is the D term for 
the SM gauge interactions i. Let us discuss the vacuum structure in this potential in the 
following. 
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First, this potential has global SUSY conserving minima at 



S = -, 1/1 = 1/1 = ^, (4) 

where the phases of / and / are constrained as arg(//) = 0. In this paper, for justifying 
the expansion of Eq. ([2]) by S/ A, we impose the following condition. 

In the case that this condition is not satisfied, the global minimum (jl]) is not determined 
within validity of this potential and dependent on UV models, so that we could not discuss 
the stability of the vacuum. 

Next, this potential also has a local minimum at the origin of the field space, 5* = 
/ = / = 0. SUSY is broken in this vacuum due to non-zero F term, Fs = /i^. The 
mass squared for the scalar component of the singlet field S is given by m| = /i^/A^. The 
stability of this metastable vacuum becomes weaker for larger A, since the potential for S 
direction becomes flat. 

The masses for the fermionic and scalar components of the messenger fields are given 
by Mmess = M and mj^ = ± k/x^, respectively, around the origin. This minimum is 
metastable when the following condition is satisfied, 

> Ki?. (6) 

When is close to k/i^, the potential for the direction along / = / is fiat and this 
metastable vacuum becomes unstable. 

The spontaneous SUSY breaking is link to the U{1)r symmetry [12]. However, this 
symmetry is approximate in this model since the SUSY mass terms of messenger fields 
violate it. The radiative correction from messenger fields destabilizes the SUSY breaking 
metastable vacuum. Thus, it is essential to include the radiative correction in discussing 
the stability of the metastable vacuum. 

The messenger fields induce the following CW effective potential. 



(\M -kS\ 



167r2 



(7) 



where 



F{x) = ^[{l + xf log(l + xf + (1 - xf log(l - xf - 2x^ log x^] . (8) 

Here, Ni is the number of messenger fields that propagate in the loop diagrams, and Ni = 5 
in our case. In the second line, the CW potential is expanded by S around the origin only 
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with the leading-order contribution about k/i^/M^, and constant terms are neglected. The 
field value of S on the local minimum is shifted toward the positive direction, and m| and 
mj^ become smaller. Thus, the CW potential makes the local minimum more unstable. 
Nevertheless, as long as A is small enough to satisfy 

M>^A, (9) 

the local minimum remains and exists near the origin: (S) ~ 5(/t^A/47rM)^M/K. At this 
false vacuum, the condition ([6]) does not significantly change, of course, as far as {S) is 
small enough. 

Now we discuss the relation between the gravitino mass and the stability of the 
metastable SUSY breaking vacuum. First we give the gravitino mass, normalized by 
the squark mass, in this model. The SUSY breaking is mediated to the visible sector by 
messenger fields. The gaugino and sfermion masses are induced at the messenger scale hm 
as 

2 



M|(,„, . pNC. {^J .pNC. (-^)^ (^)\ ,U) 

where Cj is the quadratic Casimir, A^ is the number of the messenger pair, and A^ = 1 
for our minimal case. We assumed M ^ k,{S) and Fs — /i^ in the last form. Then, the 
gravitino mass is given by 

™3..^.I0evf^Vf^V, ,12) 



VSMpi \2TeVj \KfiJ 

where Mpi = 2.4 x 10^^ GeV is the reduced Planck mass. Here the mass scale is normalized 
by the squark mass Mg calculated by Eq. ( ITTi) . 

Small M/ufi is required in Eq. (fT2l) for the ultra light gravitino. However, the small 
value of M/ k/j, may lead to the tachyonic messengers, and the large value of k also desta- 
bilizes the local minimum once radiative corrections are included. The constraints ([6]) and 
([9]) lead to the lowerbound on the value of M/ nfj, as 

f^V>^^. (13) 

Thus, the ultra light gravitino m3/2 ^ 16 eV is marginal. 

In this estimate we omitted all possible 0{1) coefficients. The precise value of the 
minimum of M/ Kfi is essential to discuss the possibility of the ultra light gravitino in this 
model. In next section we evaluate the decay rate of the false vacuum numerically and 
constrain the parameter space of this model. 
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3 Stability of metastable vacuum 



As explained in the previous section, the stabihty of the SUSY breaking false vacuum is 
weaker for smaller gravitino mass. Hence, we should carefully estimate the transition rate 
of the false vacuum to the true ones. It can be evaluated in the semiclassical techniques 
[To] . In this method, the transition rate of false vacua described by scalar fields (pi is cal- 
culated by the Euclidean path integral Jldcpi] e~'^^t<Ai]^ which is dominated by the so-called 
"bounce" configuration (pi = (j)i{x). This bounce is a stationary point of the Euclidean 
action S^i^i], that is 6SE[(l>i] = 0, and it is known to be an 0(4) symmetric solution [13] . 
Once we find the bounce solution, the transition rate per unit volume is estimated as 

T/V ~ A e~^. (14) 

The prefactor A has the fourth power of the typical scale in the potential and 

B = SEmr)]-SE[<l>{], (15) 

where 0f are constant configurations that continue to stay at the false vacuum. 

Several complex scalar fields exist in the present model. Hence, we should estimate 
Se by finding the bounce solution in a numerical method. In the following, we explain 
our calculation. For simplicity, we take all fields real and / = /. The D term potential 
for messenger fields makes this configuration |/| = |/|, and the potential for the direction 
along / = / becomes unstable when S is increased. Then, defining 0i = / = / and 
02 = S, we have the following 0(4) symmetric Euclidean action, 

drr' ^ 2 V rfr 
1=1 



(16) 



where ki = 4, k2 = 2 and V{(j)i) is the scalar potential, including the CW potential, that is 
the sum of Eqs. and (^^. Then the bounce solution 0j is given by solving the following 
equations of motion, 

k-^^L^i + k-— = —((h) (17) 
* dr'^ V dr d(hi * ' 



with boundary conditions. 



f d - 
lim (j)i{r) = (pj, = 0. 



+ 00 



Since the prefactor A is approximately given as A ~ fi'^, B > 400 is required to make the 
lifetime of the false vacuum much longer than the age of the universe. 

To solve the above equations and find the bounce configuration numerically, we used 
the calculation method discussed in Ref. [T^. We discretized Eqs. f lT7|) and ffTS]) as 



h.j+i - 2(pij + (pij-i a - I (pij+i - (pij-i _ dV 

(Ar)2 + ' (j + Aj)Ar 2Ar ~ d(Pi 



(20) 
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A//i = 5,Af//i = 2,K = 0.5 4.5 




Figure 1: A typical shape of potential (left) and dependence of bounce solution 

(f)i{r) on r (right). Parameters are taken as A//i = 5, M/fi = 2 and k = 0.5. A trajectory 
of (f)i{r) is shown by a dotted line in the left figure. Height of the potential in the figure is 
scaled by fi"^. 



where (j = 0, ■ ■ ■ , n) represents the field value of 0j at the j-th lattice site and Ar is 
the lattice size. These numerical parameters are taken to be n = 200 and Ar = 0.25//i in 
our calculation. Parameters a and Aj in Eq. (fT9l) are used to find the bounce configuration 
in this method. At first we found the undamped solution with a = 1, which can be easily 
obtained by an improved potential . Then we gradually increase a to 4 and decrease A? 



to with solving the above equations, and finally the desired bounce solution is obtaineq^ 



From this bounce configuration, the coefficient B is calculated by the following discretized 
Euclidean action, 



n-l 
j=0 



J + 



Ar^ 



E 



Ar 



+ 



(22) 



In Fig. [H we show a typical potential shape (left) and the corresponding bounce con- 
figuration derived by the above procedure (right) for A//i = 5, M/yU = 2 and k = 0.5. In 
the left figure, a trajectory between the false vacuum and true one is shown by a dotted 
line. In the right figure, the trajectories of 0i (= /) and 02 (= S) are shown by solid and 
dotted lines, respectively. They are rather simple in most of the parameter space as long 



'^'^ In addition, after one bounce solution cj)* is obtained for a parameter set, new bounce solutions for the 
new parameter sets are also iteratively derived from it. In this derivation the right-hand side of Eq. 
is expanded as 



dV dV 



2 

E( 

fe=i 



d(t>id4>k 



^i,7")j 



(21) 



to linearized the system of equations. 
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as the conditions Eqs. ([5]), ([6]) and ([9]) are satisfied. First 02 moves toward the positive 
direction and 0i is almost unchanged. Then around (p2 — {M — / k, the messenger 
field becomes tachyonic and both fields start to go toward one of the true vacua. 

Next, we discuss the parameter dependence of the decay rate of the false vacuum. The 
potential V{(j)i) are described by 4 parameters: A, M, fi, and k. Since the overall mass 
scale is irrelevant to the calculation of B, we normalized all dimensional parameters by 
/i. In Fig. [2] (left), we show the value of S as a function of M/k/i and k, for A//i = 5. In 
the lower-left region the vacuum is unstable along the /(/) direction since the condition, 
> Kfi^, in Eq. ([6]) is not satisfied. In the upper-left region the CW potential is not 
negligible and the vacuum is unstable due to the violation of the condition in Eq. ([9]). 
From this figure, it is found that B is determined mainly by the value of M/kh except 
around the boundaries given in Eqs. ([6]) and ([9]). 

In Fig. [2] (right), the allowed regions by the constraint B > 400 are shown for several 
values of A//i. The allowed regions are the right side of solid lines. The minimum allowed 
value for M/ nfi is almost determined by the intersecting point of the stability conditions, 
Eqs. (E]) and ([9]) for fixed A//i. Hence, the bound is roughly given by Eq. ( fT3|l . The differ- 
ence between Eq. ( IT3i) and the bound derived by Fig. [2] is explained by two contributions. 
First, the numerical calculation of the transition rate cuts some region near this bound 
as shown in Fig. [2] (left). The other is the approximation neglecting 0{1) factors in the 
calculation of Eq. ([9]). The rest condition ([5]) gives upperbound on M/hzfi with a constant 
A/fi. Hence, this condition is irrelevant to the lowerbound on the gravitino mass. 

The parameter dependence of B is understood from an approximate estimate based 
on Ref. |TT], where the vacuum transition rate is calculated for a triangle-shaped poten- 
tial of one real scalar field. Fortunately, the trajectory of the bounce solution 0j(r) is 
rather simple as depicted in Fig. [H except around the boundaries, Eqs. ([6]) and Qj. We 
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Figure 3: Constraints on M/nfi as functions of A//i (left). Solid line is lowerbound on 
M/nfi, while dotted lines show M/kA = 1,1/2,1/3 from left to right. Upperbound for 
the squark mass as a function of the gravitino mass (right). We take A/yU = 2, 3, 5, 10, 20. 
Allowed regions are below the lines. 



approximated the potential along this trajectory by a triangle form and derived as 

B ^ 23 



where this estimate is applicable for thin wall region^i. This approximated estimate 
explains some features of numerical results. First, the value of B depends mainly on the 
combination of M/ Kfi and the power of M/ Kfi is rather large. Second, the dependence of 
i? on A is rather weak. This implies that B does not change so much for different values 
of A as long as M/ nfi are fixed, although the allowed region are changed. 

In Fig. [3] (left), we summarized the allowed minimum values of M/k/i dependent on 
A//i. The solid line shows the lowerbound on M/k/i and the dotted lines are M/kA = 
1,1/2 and 1/3, respectively. For larger values of M/kA, the description as an effective 
low-energy theory becomes inadequate and higher order terms neglected in Eq. ([2]) become 
important. In these cases, UV models are needed for complete analysis. 

The lowerbound on M/ufi can be translated to the constraints on the gravitino and 
squark mass through Eq. f|T2l) . In Fig. [3] (right), we show the allowed region for the squark 
and gravitino masses for several values of A/yU, where the allowed regions are below lines. 
From this figure, the upperbound on the squark mass for fixed value of is obtained. 
It is found that the relatively light squark mass Mg < 1800 GeV are required to have 
^3/2 ^ 16 eV even if A//i = 2. When A is larger, the constraint is more severe. The 
ultra light gravitino in this model may be checked, since the mass range is accessible 
in the LHC experiments. This result is derived by using the approximation of the last 



Our estimate is different from that in Ref. [S], since the result in Ref. "E^ can be reahzed using a 
different assumption near the boundaries. 
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form of Eq. (fT2!) . The upperbound on the squark mass is shghtly relaxed without this 
approximation. Nevertheless, we checked that the difference is about 10-20 %. 

The LEP experiments constrain the SM-like Higgs boson mass as rrih > 114.4 GeV 
[15] . For the large enough Higgs boson mass, the stop mass should be larger than about 
1000 GeV to have a large radiative correction through top and stop loop diagrams. This 
suggests that the cutoff scale is low enough to satisfy A/yU < 20 as in Fig. [31 This is 
reduced to the bound, A < 10^ GeV, with /i ~ 200 TeV(m3/2/10 eV)^/2_ jhis gives a 
strong implication for UV models. 

Finally we briefly discuss the case of non- minimal messenger fields (A^ > 1), which 
is a quite interesting situation, since it may be possible to measure the gravitino mass 
in future collider experiments [16]. If A^ copies of 5 + 5* are introduced, the gluino and 
squark masses become larger and scale like A^ and V^, respectively. On the other hand, 
the SUSY breaking vacuum becomes more unstable, since the CW potential in Eq. ([7]) 
is proportional to A^. A factor N^^^ appears in the lowerbound on M/k/i for a fixed 
A/fi. Then the upperbounds on the gluino and squark masses scale like N^/^ and N^^^ 
respectively. Hence, squark may be still detected at the LHC experiments even for A^ > 1. 

4 Summary 

In this paper we have discussed the possibility that the ultra light gravitino with m3/2 ^ 
16 eV is realized in the GMSB model with the metastable SUSY breaking vacuum. In 
our numerical calculation of the decay rate of the false vacuum, the Coleman- Weinberg 
potential is included in the scalar potential, since it destabilizes the metastable SUSY 
breaking vacuum. We have found that when this model predicts a rather light gravitino, 
an upperbound on the squark mass Mg < 1800 GeV at most is obtained for the minimal 
messenger fields from the stability of the vacuum. They are detectable in the LHC exper- 
iments, and hence this scenario may be checked. Even for the cases of A^ > 1, squarks are 
still relatively light. We also found that to have the ultra light gravitino the cutoff scale 
must be smaller than 10'' GeV, which is an important requirement for UV models that 
realize this scenario. 
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